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Abstract 

We develop a wavepacket approach to the diffraction of charged particles by a thin material target 
and we use the de Broglie-Bohm quantum trajectories to study various phenomena in this context. 
We construct a particle wave function model given as the sum of two terms t/r - iff mgoing+^ outgoing, 
each having a wavepacket form with longitudinal and transverse quantum coherence lengths both 
finite. We find the form of the separator, i.e. the limit between the domains of prevalence of the 
ingoing and outgoing quantum flow. The structure of the quantum-mechanical currents in the 
neighborhood of the separator implies the formation of an array of quantum vortices (nodal point 
- X point complexes). The X point gives rise to stable and unstable manifolds, whose directions 
determine the scattering of the de Broglie - Bohm trajectories. We show how the deformation 
of the separatior near Bragg angles explains the emergence of a diffraction pattern by the de 
Broglie - Bohm trajectories. We calculate the arrival time distributions for particles scattered 
at different angles. A main prediction is that the arrival time distributions have a dispersion 
proportional to v^'x the largest of the longitudinal and transverse coherence lengths, where vq 
is the mean velocity of incident particles. We also calculate time-of-flight differences AT" for 
particles scattered in different angles. The predictions of the de Broglie - Bohm theory for AT" 
turn to be different from estimates of the same quantity using other theories on time observables 
like the sum-over-histories or the Kijowski approach. We propose an experimental setup aiming 
to test such predictions. Finally, we explore the semiclassical limit of short wavelength and short 
quantum coherence lengths, and demonstrate how, in this case, results with the de Broglie - 
Bohm trajectories are similar to the classical results of Rutherford scattering. 
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1. Introduction 

The de Broglie - Bohm quantum trajectories [1] [2] [3] [4] have been considered as an in- 
terpretational tool in a number of recent applications (see [5] [6] [7] for reviews), since they can 
offer new insight into a variety of complex quantum phenomena. According to the de Broglie- 
Bohm theory, to any wavefunction i/'(ri, r2, . . . , Vj^i, t) describing a A^-particle system, we can 
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associate a set of 'quantum trajectories'. One trajectory is defined by the initial conditions 
(ri(0), r2(0), . . . , riv(O)) and by the 'pilot wave' equations of motion 



-j- = -Im{-^\ i = l,...N (1) 
at rrii ij/ 

where m, are the particle masses and h is Planck's constant. The equations of motion (1) im- 
ply the continuity equation for the probability density p(ri, r2, .. .,YN,t) - \4'{r\ , r2, . . . , rj/, 
In particular, in a one-particle system we can choose many different initial conditions corre- 
sponding to an initial density p(x, 0) - \4i{r, 0)p. Then, the pilot-wave equations guarantee the 
preservation of Horn's rule p(r, f) - |i/'(r, f)P at all subsequent times t. Furthermore, the de 
Broglie - Bohm trajectories are equivalent to the stream lines of the quantum probability current 
j = (?i/2m;)(i/''V(/' - if/Vi//*). Thus, the Bohmian approach yields practically equivalent results to 
Madelung's quantum hydrodynamics [8]. 

The de Broglie - Bohm theory has been discussed extensively from the point of view of its 
relevance as a consistent interpretation of quantum mechanics (e.g. [4], [9]; see [10] for an 
extended list of references). However, the employment of the de Broglie - Bohm trajectories 
has been proven useful also in many practical aspects of the study of quantum systems. Some 
modern applications are: 

i) Visualization of quantum processes: examples are barrier penetration or the quantum tun- 
neUng effect [11] [12] [13] [14], the (particle) two-slit experiment [15], balUstic transport through 
'quantum wires' [16][17], molecular dynamics [18], dynamics in nonlinear systems with classi- 
cal focal points or caustics [19], and rotational or atom-surface scattering [20] [21][22]. 

ii) Lagrangian solvers of Schrodinger equation via swarms of evolving Bohmian trajectories 
(see [5] for a comprehensive review, as well as [23] [21] [22] [24]). The interest in this method 
lies in that, instead of solving Schrodinger's equation first, one uses a step-by-step procedure to 
calculate the trajectories via Newton's second order equations of motion in a potential 

Uir,t) = Vir,t) + Qir,t) (2) 

where Q{r, t) is the 'quantum potential', caused by the wavefunction ij/: 

Using the information of the initial value of the wavefunction as well as the evolution of the 
quantum trajectories, the wavefunction can then be determined at any subsequent time step. 

in) Dynamical origin of the quantum relaxation [25] [26] [27] [28] [29]. The de Broglie - 
Bohm theory offers a justification of Bom's rule p = \\lr^, since it predicts that, under some con- 
ditions, the quantum trajectories lead to an asymptotic (in time) approach towards this rule even 
if it was initially allowed that Pmiuai + initial?' ■ It should be noted that not all choices of pi„itiai 
are guaranteed to lead to quantum relaxation, and counter-examples can be found, for reasons 
explained in [26]. The arguments used in that paper to explain the suppression of the quantum 
relaxation effect in the two-sht experiment apply also to many other cases (see e.g. [30]). In par- 
ticular, a necessary condition for quantum relaxation to take place is that the trajectories should 
exhibit chaotic behavior (see [25] [26]); however, even this condition is not sufficient (see [31]). 
The problem of chaos in the de Broglie - Bohm theory has been studied extensively (indicative 
references are [32] [33] [34] [35] [3 6] [3 7] [3 8] [3 9] [40][41][42][43] [44][45][46][47]). We have 
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worked on this problem in [26] [48] [49] [50] [51]. Our main result was that chaos is due to the 
presence of moving quantum vortices forming 'nodal point - X-point complexes' ([48] [49] [50]; 
see also [45][46]). Quantitative studies of chaos and of the effects of vortices are presented in 
[38] [47] [50]. In particular, in [50] we made a theoretical analysis of the dependence of Lya- 
punov exponents of the quantum trajectories on the size and speed of the quantum vortices, thus 
explaining numerical results found in [48] and [49]. Furthermore, in [26] we gave examples of 
systems which do or do not exhibit quantum relaxation, depending on whether or not their under- 
lying trajectories are chaotic. It should be emphasized that, besides chaos, the quantum vortices 
play a key role in a variety of quantum dynamical phenomena (e.g. [52][1 1] [53][21][22]). 

iv) Arrival times and times of flight. In the traditional formulation of quantum mechanics 
time is only a parameter in Schrodinger's equation, since by a theorem of Pauli [54] no defi- 
nition of a self-adjoint time-operator consistent with all axioms of quantum mechanics can be 
given in a system with energy spectrum bounded from below. Time, however, is an experimental 
observable. Various approaches in the literature, reviewed in [55] [56], have addressed the ques- 
tion of a consistent definition of quantum probability distributions for time observables. Besides 
the Bohmian approach, two other approaches are: a) the 'sum-over-histories' approach [57] [58] 
based on Feynman paths, and b) the approach of Kijowski [59], based on the definition of quan- 
tum states acted upon by the so-called 'Bohm-Aharonov operator' (see [55]). On the other hand, 
the de Broglie - Bohm approach gives a straightforward answer to this problem, since the time 
needed to connect any two points along a quantum trajectory is a well defined quantity (see [60] 
[61]). 

Regarding this latter point, a key remark that will concern us in the sequel is that a consistent 
definition of the arrival times, that would allow in principle for a comparison of the various 
approaches in specific quantum systems, is only possible provided that the initial wavefunction 
is localized in space, i.e. it is described by a wavepacket model. 

Being motivated by the latter remark, in the present paper we present a study of the de Broglie 
- Bohm trajectories in a wavepacket model referring to a quantum phenomenon that has played 
a fundamental role in the development of quantum mechanics, namely the diffraction of charged 
particles (e.g. electrons or ions) by a thin material target. 

A theoretical study on the quantum scattering problem in the framework of the de Broghe - 
Bohm approach has been presented in the series of works [62] [63] [64] [65] [66]. These studies 
refer to the establishment of the rules of scattering probabilities using the 'flux across surfaces' 
theorem adapted to the concept of quantum trajectories. However, they do not deal with the form 
of the quantum trajectories or the emergence of difiraction patterns under specific scattering po- 
tentials. A numerical simulation of Rutherford scattering by a single nucleus has been presented 
in [67], while in [23] the phenomena of atom-surface scattering as well as neutron difiraction by 
slits are considered, which share some common features, but also important diff'erences, with our 
problem. 

In the present paper we make a detailed study of the de Broglie - Bohm trajectories in the 
context of a wavepacket model of charged particle diffraction, by first investigating the form of 
the quantum currents corresponding to various cases of this model. These cases are diversified 
one from the other by the different quantitative relations characterizing the so-called quantum 
coherence lengths in the longitudinal and transverse directions of the charged particle beam. 
This is necessary in order to be able to compare the results corresponding to possibly different 
experimental realizations of a charged particle beam, as e.g. in the case of electrons produced 
either by thermionic or by a cold-field emission processes. 

Our present study completes in a substantial way the study initiated in a previous paper 
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Figure 1: The basic setup of the problem under study. A source (S) emits charged particles described by an 'ingoing' 
wavefunction having the form of a wavepacket with dispersions / in the longitudinal direction (z-axis = direction of 
incidence to a thin material target (C) placed at the center O of the coordinate system), and D in the transverse direction. 
After scattering, some particles arrive at detectors D, placed at equal distances from O and various angles 8j. The 
wavefunction is assumed to have axial symmetry (around the z-axis), thus the figure corresponds to any meridian plane. 
Various other symbols are explained in the text. 

of ours [51], in which we implemented the de Broglie - Bohm approach in the case of electron 
diffraction through a thin crystal. In that study, however, we assumed a planar wave model for the 
propagation of the electron wavefunction in the longitudinal direction. In contrast, in the present 
paper we assume instead a finite longitudinal quantum coherence length. This assumption leads 
to a number of crucial new elements with respect to [5 1 ] . In fact, in order to achieve our goal 
we derive a wavefunction model by a refined implementation of basic scattering theory, so as 
to account for a fully-localized in space description of scattering. The derivation of this model 
presents its own interest, and it is exposed in detail in section 2. 

The structure of the paper is as follows: after the derivation of the basic wavefunction model 
in section 2, we pass to a study of the quantum trajectories in section 3. Here the emphasis 
is on the influence upon the trajectories of quantum vortices, whose appearance and role in 
this problem are explicitly discussed. In fact, we show that the quantum vortices appear in 
the transition zone from a domain of predominance of the ingoing wavefunction to a domain 
of predominance of the outgoing wavefunction. Inside this zone we can define a locus called 
separator, which plays a key role in the interpretation of the scattering process via the quantum 
trajectories. In section 4 we study the arrival times of diffracted particles to detectors placed in 
various scattering angles. A main outcome of this study is that it is possible to propose a feasible 
experimental test probing the predictions of the Bohmian theory about the particles' arrival times. 
In section 5 we discuss separately the 'semi-classical' case of particles with a large mass and and 
a small de Broglie wavelength, applicable e.g. to a-particle or ion scattering, since this case 
exhibits some special features in comparison to the case of electron diffraction. Finally, section 
6 summarizes the main conclusions of the present study. 
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2. Modelling of the wavefunction 



We consider a cylindrical beam of particles of mass m and charge Ziq^ incident on a thin 
material target. We set the center of the target as the origin of our coordinate system of reference, 
and use both cylindrical coordinates (z, R, 0) and spherical coordinates (r, 6, (/>). The z-axis is the 
beam's main axis, R denotes cylindrical radius transversally to z, <f> is the azimuth, r = (z^ +R^y^^ 
and = tan '(/J/z) (see Figure 1, schematic). 

A basic form of diffraction theory for charged particles, reviewed e.g. in [68], assumes that 
the incident waves are planar. As explained in the introduction, here instead we are interested in 
a wavepacket approach. Focusing only on elastic scattering phenomena, the latter approach can 
be obtained by a refinement of the basic theory as follows: 

The potential felt by a charged particle approaching the target can be considered as the sum 
of the individual potential terms generated by every atom in the target: 

N 

V(r) = 2£^(r-r,) . (4) 

7=1 

where denotes the position of j-th atom in the lattice of the target (this position exhibits some 
statistical fluctuations due to thermal oscillations etc; the effect of these fluctuations is discussed 
later in this section). As a model for the function U, we can adopt a screened Coulomb potential 

1 ZiZfl? exp(-|r - ril/rn) 

f/(r - rj) = , , ' (5) 

4neo |r - rjl 

(eo = vacuum dielectric constant), where Z is the nuclear charge, and ro is a constant representing 
a charge screening range within the atoms, whose value is of the order of the atomic size. 

Particles being scattered by the target can be described by a wavefunction given as a super- 
position of eigenfunctions 



(6) 



where c(k) are Fourier coefficients, and ^k(r) are scattering eigenfunctions, i.e. solutions of the 
time-independent Schrodinger's equation 

_^v20 + y(r)0 = £0 (7) 

Zm 

with V chosen as in (4) and £ > 0. The difi'erent solutions ^ = (pkare labeled by their wavevectors 
k of modulus A: =| k |= {2mE)^l^ 1%, where E > is the energy associated with one eigenstate. 
Bom's approximation can be used to obtain an approximative formula for 0k- We thus write 

0k = 0O,k + <;*l,k + (plM- + • • • (8) 

where 0o,k = e'^ = 0(l) is the solution of Eq.(7) for the free particle problem (y(r) = 0), 
while 0i_k = 0{V), (/>2,k = 0(V^) etc (assuming that V small compared to the particles' energies). 
The above series are meaningful at aU points of space excluding a set of balls of radius a few 
times ro around every one of the atoms in the target. Spherical harmonic expansions (see e.g. 
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[69]) provide a more accurate representation of the solution inside such balls, but their use is 
cumbersome while practically unnecessary in the context of the present study. 

A step by step determination of the series terms in (8) can be obtained via the recursive 
formula 



-ir-^ <f>nM + V<pn-l,k = E(j)n,Vi 

2m 



2m 



■<Pn.k 



(9) 



All essential phenomena discussed below are present already in the solutions including just the 
two first terms (pk - <Po,k + <Pi,k- From Eq.(9) for n = 1 we find: 



<;*i,k(r) 



m r 
"2^ Ja] 



*|r-r'| 



all space 



Ir - r'l 



N 

Jk-V V 



1 ZiZ^2g-|r'-rj|/ro 



Ir' - nl 



(10) 



The integral in (10) can be estimated using standard approximations of scattering theory. We 
then find 



0k(r) 



Jk-r 



ZiZq^ m 



( N 



gik\T-T j\gik-rj 



■rj\i2k^sin\Aej/2) + l/2rl)) 



(11) 



where A6j denotes the angle between the vectors k and r - rj. 
Substituting Eq.(l 1) into Eq.(6) we have 



(23 

ZiZql m 



(12) 



g*|r-rjlg!k-rj 



[j^,\r-rj \(2k^ sm\Aej/2) + l/2rl) ) 



The problem of defining il/{v, t) is now restricted to making an appropriate choice for the co- 
efficients c(k). The latter are determined by the Fourier transform of the initial wavefunction 
il/{r, f = 0). In the wavepacket approach, the initial wavefunction is localized around the source, 
i.e. far from the target. Hence we can set ^(r, ; = 0) ^ >/fingomgi^< t = 0), where ifrmgoingi^, t = 0) 
represents a wavepacket moving in the z-direction towards the target with some velocity vq. A 
Gaussian wavepacket of this form corresponds (in momentum space) to the choice 



c(k) 



1 



1 



exp 



2al 



2 \ 



exp 



(13) 



In (13), (kx,ky,ky) are the Cartesian components of k, zo - -Iq is the initial position of the 
center of the wavepacket along the z-axis, and = mvo/h. The quantities cry, crj_ are the 
longitudinal and transverse dispersions of the wavepacket in momentum space. These correspond 
to dispersions in position space given by Z = cr'^^ and D = cr^^ The quantities I and D are 
hereafter called the longitudinal and transverse quantum coherence length respectively. Eq.(12) 
now takes the form 

ij/ir, t) = iJ/i„goi„gir, t) + )J/outgoingir, t) (14) 



where 



(z + k- ^tf 



with 



Bit) 



1 



D 



I 



1/2 



exp 



ikolo ■ 



2m ' 



„3/4 ^1)2 + ihtim j\P + iht/m j 

The function (jjomgomg corresponds to the second integral in (12). An expHcit expression for 
this function can only be found by adopting some further approximations. First, we consider 
fast-moving wavepackets, for which » niax(crj_,cr||) as well as ko » l/''o- Then, in the 
denominator of the second integrand in (12): i) the term l/2rQ can be ignored, and ii) we use 
the approximation Ijk^ ^ Second, at all distances r » ry we have that the angles A^y are 
approximately equal one to the other and to the angle 9 (which is equal to the angle between the 
vectors r and ko = (0, 0, fco)- Finally, we set | r - r j |a; r - rj ■ n + rj/(2r) in the exponential 
argument of (12), where n = (sin cos ^, sin sin cos 0) (this is necessary in order to retain 
all terms whose phase has a substantially non-zero value), while we set | r - rj |« r in the 
denominator of the integrands in (12). Using these approximations, we find 

ZiZql m 1 expiikor) 



outgoing 



Bit)- 



z 

7=1 



4;reo 2kl?.m^{ei2) r 
exp (jA;o[-rj -11-1-2^ + rj/(2r)]) 



(16) 



exp 



(r + /o - vo? - Tj 



n + ^7 + 2r 



1)2 



2(£)2 + m.-) 



2(/^ + f) 



where vq = hko/m represents the mean velocity of a particle with wavenumber^o- 

At distances r closer to the target than the maximum of the two coherence lengths D,l, 
the prefactor f{r,6) = l/{2kQsm^i6/2)r) in Eq.(16) is no longer accurate. In order to be able 
to perform some numerical calculations of de Broglie - Bohm trajectories, after numerically 
simulating the sums appearing in (12) we found by trial a fitting model that represents reasonably 
well the modifications of /(r, 6) close to the target. This reads: 

fir, 0) = ko^ [c^D sin 6 + (cjo^ sin^ 6 + - 2rc4D sin 6 + cjo^)^'^ - r cos 0^ (17) 

where C3 and C4 are fitting constants determined by comparison of Eq.(17) to the results of the 
numerical simulation of f(r,ff) (in all simulations below we set C3 = 0.3, C4 = 0.8). Il is to be 
stressed that Eq.(17) correctly recovers the asymptotic form / 
large. 

The outgoing wavefunction now takes the form 

B(t)ZiZql}77 



\l{2kl sm^{ei2)r) when r is 



4^, 



outgoing 



e'''>'y{r,0)S,ff(ko,r,t) 



(18) 



where the quantity Se/fiko; r) is called hereafter the 'eff"ective Fraunhofi'er function' (in analogy 
with the 'far field' diflraction limit in wave optics, see [70]). This is given by 



N 

Seff(ko,r,t) = Yj 



exp [ikoi-Tj ■n + Zj + rj/(2r))) 



(19) 



The physical significance of the function Sef/iko; r, t) is that it sums the contributions of all the 
atoms in the target which act as sources of partial outgoing waves, whose superposition forms 
outgoing- Furthermore, the function Se/f accounts for the formation of a diffraction pattern, 
which, for given n, r, arises by the coherent contributions of all atoms in the target whose phasors 
exp[i7jo(-rj -11 + Zj + rj/i2r))] are nearly parallel one to the other. As in [51], we consider the 
simplest example of a cubic lattice structure of the target 

Vj = {Hx, fly, n^a + AaUjit), 

(20) 

e (__,_) x (--,-) x (--, -) 

where i) a is the lattice constant (equal to the length of one side of the primitive cell) ii) Aa is the 
amplitude of some random oscillations (due to thermal or recoil motions; Aa is taken equal to a 
small fraction of a) and uj = (uj^^, My , , My ,) are random variables with a uniform distribution in the 
intervals [-0.5, 0.5] (the random oscillations introduce a so-called Debye- Waller effect, analyzed 
in [51]; here, for simplicity, we ignore modifications on the wavefunction due to this effect), iii) 
The number of atoms A^, in the z-direction is A^, - d/a, where d is the target thickness, and iv) 
the value of N± is of order A'j, - OiD/a), due to the Gaussian factor exp(-/?J/(2Z)^ + iht/m)) 
in Eq.(19) which can be approximated by « 1 for all In^l < Nj_/2 and \ny\ < N±/2, and by for 
> A^±/2 or |«vl > A^±/2 (for typical magnitudes of D the inequality » ht/m holds for all 
times t of interest in our study, see below). We now distinguish the following cases: 

2.1. l» D» a 

When the longitudinal coherence length I is larger than the transverse coherence length D, a 
simple modeling of the sum in Eq.(19) becomes possible at all distances r > D. Ignoring first 

the random fluctuations in (20) (i.e. setting Aa = 0), we have r» | - rj ■n + Zj + ^\ whereby it 
follows that 

S effiko, r, t) « exp - — T — 

\ 2(P + iht/m) 

iko[-anx sin 6 cos - any sin sin H ] | (21) 



n,=-NJ2ny=-NJ2 

^ exp I jIo [( 1 - COS 0)n^a + -|— ] j 



n,=-NJ2 

For a random choice of ko, 6, (p, the total number of contributing atoms in the sums of Eq.(21) is 
of the order A'^ ~ Nj_N^ = D^d/a^. Furthermore, the A' phasors have an effectively random phase. 
Thus, the total sum is of order A'^^''^, and we are lead to the simple estimate 

Seff ~ Dd^l^la^l^ exp[-(r + /q - voOV(2(/^ + jfe/m))] 

called, hereafter, the 'diffuse term' of the effective Fraunhofer function. Using this term we have 



B{t)ZiZqlmDd'^l^p^l^ 
8 



^ — \-/ — 1 — ~ie r- 

Y outgoing ~ ^ exp 



(r + /()- VQf)^ ' 
2{P + ^) 



f{r,ey'^' (22) 



where p = is the number density of the atoms in the target. 

The model (22) is used in a number of numerical simulations below. It describes a radial 
pulse propagating outwards with speed vq, which emerges from the center at the time /q/vq. It 
should be stressed, however, that Eq.(22) requires modifications close to particular angles where 
the phasors in the sums of (21) are added coherently. This will be examined in subsection 3.3. 
Further modifications are required when / and D become comparable or smaller than the inter- 
atomic distance a in the target. This case is examined in section 5. 

2.2. D»l»a 

The modeling of S ej-j is a more subtle problem if the transverse coherence length is larger 
than the longitudinal coherence length. Without loss of generality, we can consider a fixed merid- 
ian plane defined e.g. by the angle = 0. We set ^ = r + Iq - vot and we make the change of 
variables u = -xsinO + (x^ + y^)/(2r), R - (x^ + y^Y^^. Considering now a random choice of 
^0, 0, and ignoring the quantity iHt/m, for a given value of ^ the sum over the variables rix, riy 
in Eq.(19) can by approximated by a sum over a domain of values of (nx, riy) such that u{nx, riy) 
belongs to a ball of radius I around f . The total number of contributing atoms is then of order 
N ~ fid/a^. If i) we consider the sum over the phasors as a sum of random numbers (yielding a 
total magnitude ~ A'^'^^), and ii) we substitute the second exponential in (19) by a delta function 
around ^, we find 5^// « {ld^^^/a^^^)I, where 



Jiu, Kje-^'i^^'diu + OdudR (23) 

and J{u,R) is the determinant of the Jacobian matrix of the transformation {x,y) — > {u,R). We 
find 



■JR (24) 



where 7?m,„ - \r(sm0 - -^sin^ 6 - 2^/r)|, R^^ - r(sm0 + -^sin^ - 2^/r). Eq.(24) yields non- 
zero values of the integral / below a cut-off radius r < (vq/ - /o)/(l - 5 sin^ 6*). However, the 
asymptotic behavior of / when r is large is found by noticing that « ^/sin6 and Rmax — > °° 
in this limit. We then find 



The essential point to retain is that the profile of the radial outgoing pulse in the case D » / is 
a Gaussian whose dispersion is of order D. Thus, the conclusion is that, in both cases I > D or 
D > I, the outgoing wavefunction has always the form of a packet with dispersion cr^ of the same 
order as the largest of the two quantum coherence lengths, i.e. (Xr ~ max(/,Z)). 

Furthermore, we notice that in the case D » / the dispersion cr^ depends also on 0. A 
detailed investigation of the Bohmian trajectories in this case is, however, not possible from a 
numerical point of view, because the asymptotic formulae (24) and (25) are not valid at distances 
r < D, where the scattering efi'ects take place. Thus in the sequel we limit ourselves to a detailed 
investigation of the Bohmian trajectories in the case / » D, while a qualitative discussion of 
the case D » I will be made in section 4, referring to the issue of the particles' arrival time 
distribution. 
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3. Quantum trajectories 



We now discuss the main features of the de Broghe - Bohm quantum trajectories focusing on 
the case / » D. 

3.1. Separator and quantum vortices 

The form of the trajectories can be found by carefully examining the structure of the quantum 
currents j - {h/2mi){i(/*'Vil/ - i/'Vi/'*). The main remark is that, due to Eqs.(15) and (22), the 
ingoing wavefunction term (which has a Gaussian form both in the R and z directions) has a 
falUng exponential profile at large distances from the center of the Gaussian, while the outgoing 
wavefunction has a more complex form falling asymptotically as a power-law 1/r due to the 
factor /. Thus, there is an inner domain of the quantum flow where (//ingoing prevails, and an outer 
domain where ifroutgoing prevails. We call separator the boundary dehmiting the two domains. 
Formally, the separator is defined as the (time-evolving) geometric locus where 

\^ingoing\ — \^outgoitig\ (26) 

In the case where the outgoing wavefunction is given by Eq.(22), the condition (26) takes the 
form: 

where we use the approximations + iht/m ^ and F + iht/m ^ F. We note that these 
approximations hold within a range of parameter values relevant to concrete experimental setups. 
For example, assuming that incident particles have velocities of order vq ~ lO^m/s, the time 
required to travel a distance of order 10"^ - 10"^m (which is the typical size of an experimental 
setup) is of the order of / = 10 '° - lO^^s. On the other hand, the typical coherence lengths in 
experiments e.g. for electrons are of order Ijum or larger. Hence, ht/m is much smaller than 
or/2. 

The time evolution of the separator in the plane (R, z) depends now on the time evolution 
of the relative amplitude of the ingoing compared to the outgoing wave at any point of the con- 
figuration space. We note first that according to Eq.(22), the outgoing wave corresponds to a 
wavepacket with dispersion / which emerges from the center in the time interval to < t < f^, 
with to = Ho - l)/vo, = (/o + /)/vo, which is the interval during which the support of the ingo- 
ing wavepacket (moving from left to right in Figs.l and 2) essentially overlaps with the spatial 
domain occupied by the atoms in the target (see [69] for an introductory description of this phe- 
nomenon in a simple Rutherford scattering case). As indicated by Eq.(22), after its emergence 
the outgoing wave moves in all radial directions maintaining essentially its Gaussian profile, 
while its overall amphtude drops like r ' . As the outgoing wave moves outwards, it first encoun- 
ters the ingoing wavepacket at times close to to. In Fig.2a, this encounter results in a gradual 
approach of the separator towards the z-axis (the indicated times are ti - 0,t2 = 3/o/(5vo) < to, 
to < h = 6/o/(5v()) < fy). As, however, the ingoing packet moves from left to right in Fig.l, its 
center crosses the target at the time t = /q/vq. Afterwards, the ingoing wave emerges from the 
right side of the target, and its support Ues nearly completely in the semi-plane z > 0. At a still 
longer time (?4 = 9/o/(5vo)), the center of the outgoing wavepacket has traveled a distance ^ 2.51 
apart, and there is no longer any overlapping between the ingoing and outgoing wavepackets. As 
observed in Fig. 2a, a transition takes place at some time between t^ and t4, such that, before the 
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Figure 2: (a) The form of the separator at four different time snapshots t\ = 0, (2 = 3/o/(5vo), h = 6/0 /(5vo) and (4 = 
9/o/(5vo) in the model where ^ingoing is given by Eq.(15), if/outgoing is given by Eq.(22), and the parameters are Zi = -1, 
m = nie, ko = 8.877 X lO^nm"' (corresponding to electrons with energy E = 30KeV, or wavelength io = 7 x 10~'nm), 
D = lOOOnm, / = lOOOOnm (corresponding to transverse and longitudinal quantum coherence lengths 1/jm and lO/jm 
respectively), Iq = 31, Z = 79 (gold), d = 420nm, a = 0.257nm. (b) The form of the quantum current flow at the snapshot 
t=t2. 
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Figure 3: Local form of the quantum flow at the 'nodal point - X-point complex' (quantum vortex) around the nodal point 
(N) with coordinates R = 1934.42nm, z = 137.178imi in the model with parameters as in Fig.2 at the time t = /q/vq. The 
thick solid curves show the unstable (U,U') and stable (S,S') asymptotic manifolds of the X-point (X) formed imder the 
instantaneous portrait of the quantum flow. 
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transition, the separator is formed by a a pair of open curves on either side of the axis R - 0, 
while after the transition there is only one closed curve intersecting twice the axis R = Q both for 
z > and z < 0. Taking into account the cylindrical symmetry around the z-axis, the form of 
the separator in space before the transition is a cylindrical-like surface of rotation, while after the 
transition it becomes a prolate spheroidal-like surface. In fact, this time-changing surface marks 
a sharp Umit between the domains of prevalence of the axial ingoing flow and the radial outgoing 
flow, as shown in Fig.2b for the time t - t2- 

If Eq.(27) is supplemented by an equation for the phases of the ingoing and outgoing waves, 
which, for Zi < takes the form 

IcqR tan(0/2) -n = 2qn q&Z , (28) 

a simultaneous solution of Eqs.(27) and (28) defines the set of all points of the configuration 
space where the total wavefunction (Eq.(14)) becomes equal to zero. Such points are caUed 
'nodal points' . 

Around the nodal points, the quantum flow forms quantum vortices (Figure 3). The local 
form of the quantum currents in a vortex domain is very different from the general flow shown in 
Fig. 2. If we 'freeze' the time t, the instantaneous pattern formed by the vector field of quantum 
probability current j corresponds to a characteristic structure called quantum vortex, or nodal 
point - X-point complex [48] [49] [50]. That is, close to a nodal point we find a second critical 
point of the flow, where one has j = 0. This is called an 'X-point', since it can be shown that 
it is always simply unstable, i.e. there are two real eigenvalues of the matrix of the linearized 
flow around X, which are one positive and one negative. Accordingly, there are two opposite 
branches of unstable (U,U') and stable (S,S') manifolds emanating from X. On the other hand, 
the nodal point can be an attractor, center, or repellor. This determines the local form of the 
invariant manifolds U and S. It has been established theoretically [48] that, except for a set of 
very small measure, most quantum trajectories avoid the nodal point, being instead scattered 
along the asymptotic directions of the manifolds of the X-point, leading to large distances from 
the nodal point - X-point complex. Furthermore, while, in general, the motion of the nodal point 
- X-point complexes introduces chaos ([38][45] [48][50]), in the present problem this effect is 
negligible because i) the speed of vortices is extremely small (of order ~ h/(komD^) « vq), 
and ii) the quantum trajectories exhibit only a small number of encounters with nodal point - 
X-point complexes, as will be shown with numerical examples below. In conclusion, the effect 
of the nodal point - X-point complexes on the trajectories can be described as a scattering process 
without recurrences. 

For the model parameters used in Fig. 3, the size of the quantum vortex, estimated by the 
distance Rx from the nodal point to the X-point, is of the order of 10"'^m. The size of vortices in 
the present model is in fact time dependent. However, in the time interval t2 < t <t3 when there 
is essential overlapping of the ingoing and outgoing wavefunction terms, Rx is approximately 
constant and it is given by the same estimate as in the second of the equations (25) of ref.[5 1], 
namely 

1 



Rx = 



Dkl; 



(29) 



The above estimate is obtained by expanding the wavefunction around a nodal point up to terms 
of second degree in (R-Rq) and (z-Zo), where (Rq, Zo) are the coordinates of the nodal point, and 
by applying general formulae derived in [50] regarding the dependence of Rx on the coefficients 
of this local expansion. 
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3.2. Trajectories 




Figure 4: (a) A swarm of Bohmian trajectories in tlie same wavefunction model as in Fig. 2. (b) Tlie resulting radial 
distribution Pradial(>",Oi) for sixteen different angles (9, = 5° + (i - 1)10°, i = 1,2, ...,16. The near coincidence of 
all distributions Pradiali'"^ Si) and with the theoretical profile corresponding to the outgoing wavefunction model (22) 
indicates the degree of preservation of the continuity equation by the numerical trajectories of (a). 

The deflection of an orbit happens at the crossing of the separator, and it is due to the orbit 
necessarily following the flow imposed by the asymptotic manifolds of the X-points that exist 
along the separator. Figure 4a shows a swarm of 625 quantum trajectories in the model with 
same parameters as in Fig.2. The initial conditions are taken on a regular grid 25 x 25 with 
-/o - 2/ < z < -lo + 21 (where k = 31) and Z)/100 <z<4D. 

An important test of the correctness of the numerical calculations is by checking whether 
the probabilities associated with the chosen initial conditions of the quantum trajectories re- 
spect the continuity equation. To this end, setting i// ^ ^ingoing at f = 0, a volume of ini- 
tial conditions AVq - InRoAR^Azo centered around the point (zo,^o) has an associated prob- 
ability AP =1 fingomg{zQ,RQ,t - 0) p AVq- Let(zo,^o) (r, 0) be the mapping from ini- 
tial conditions to a trajectory' s coordinates at f = 21qIvq. We want to estimate the mapping 
of the probabilities Af from the volume Vq to the image of this volume under the mapping 
{zq,Ro) (r, 6). This is obtained numerically, by quadratically interpolating first the functions 
r(zo,/?o) and 0(zo,Ro) from the data available by the integration of the orbits with initial condi- 
tions on the grid of points described in the previous paragraph. The quadratic interpolation aUows 
to obtain local approximations to the functions r{zo,Ro) and 6{zo,Rq) by formulae of the form 
r = Ao+Ai(z-zo)+A2(^-^o)+A3(z-zo)(^-^o), = Bo+Bi(z-zo)+B2(R-Ro)+B3(z-zo)(R-Ro), 
where the coeflicients A,, 5; change values at every grid point. These formulae, in turn, allow to 
numerically compute the Jacobian determinant J(r, 0; zo,Ro) - ArAO/ARoAzo. FinaUy, we com- 
pute the probability function Pmdiaiir, 6) - N{9)\\j/i„{zQ, Rq, t = 0)\^RoJ{r, 6; zo, ^o). where zo, Ro 
are functions of (r, 6) and ATis a normahzation constant. 

Figure 4b shows Pradiai as function of r for 16 difl'erent values of in the interval 5° < 6 < 
165°. The fact that all curves nearly coincide implies that the numerically computed Bohmian 
trajectories respect the continuity equation of the quantum flow. 

In the analysis of the arrival times or the times of flight in section 4, use is made of the 
following information: we seek to determine, as a function of the scattering angle 6, the locus of 
aU initial conditions on the (zo.^o) plane, whereby the trajectories are eventually scattered close 
to the angle 6. 
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Figure 5: The nearly straight gray zones in the bottom left side correspond to the loci (1 or 2) of initial conditions for 
which the corresponding de Broglie - Bohm trajectories end in the angular sectors (1) = 54° ± 5° and (2) S = 134° ± 5° 
at the end of the numerical integration. The black lines on top of the gray zones (1) and (2) correspond to the fitting by 



The determination of these loci follows by approximating all the quantum trajectories as 
piecewise straight lines. Namely, from Fig. 4, it is evident that any trajectory can be considered as 
a nearly perfect horizontal line up to the point (in space and time) where the trajectory encounters 
the separator. In fact, as we can see in Fig.2, the general motion of the separator itself, as 
t increases, is downwards. That is, if we consider a ray from the center outwards along any 
fixed value of the angle 6, the separator intersects this ray at a continually decreasing value 
of r, denoted by rs(t; 6). Since all trajectories are horizontal before the encounter, we have 
that R{t) - R(()) for the z-coordinate of a trajectory with initial conditions (z({)),R(())). Then, the 
encounter takes place at the time / = tcou when R{tcou) = ^(O) = rs{tcoU', 0) sin 6. The last condition 
determines the time tcoii, which is given by 



mine R(0) as a function of Z(o)- In the limit fig^(0)/(4D^) » 1, we find an approximative formula 
by replacing the r.h.s. of Eq.(3 1) by a constant average value, i.e. 



Eq.(32). 



tcoU = 



-Z(o) +R(0)cote 



(30) 



Substituting Eq.(30) in the separator equation (27), with R = /?(0) we find: 




4D^ 



-(Z(0) + k) 
PgiO) 



R(P) = Rc + 



(32) 



where we take Rc equal to the root for R^o) of Eq.(31) when Z(o) = Zc = -lo- Eq.(32) is an 
analytical expression which gives the locus of initial conditions of trajectories that are scattered 
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close to the angle 0. Figure 5 shows lines of this form for two angles 0i - 54 , 02 - 134 along 
with the loci, at (f - 2/o/vo), formed by the final points of the trajectories scattered in a bin 
around these two angles. 



3.3. Emergence of the diffraction pattern 

As mentioned in subsection 2.1, Eq.(22) provides an approximation to the outgoing wave- 
function for nearly all sets of values (^o, (/>) except very close to combinations resulting in the 
appearance of a diffraction pattern. For specific values of this pattern can be non-axisymmetric. 
Here, however, we examine for simpUcity only the appearance of Bragg angles for which the re- 
sulting diffraction pattern is axisymmetric. This implies considering the double sum over W i, n,. 
in Eq.(21) as a sum of random phasors, while allowing for a coherent addition of the phasors in 
the second sum of (21). Eq.(19) takes the form 

Seffiko,r,t)^iD/a)e 2j e"^^^^"'"'^"'"^"^'' '^^'^K (33) 

n,=-NJ2 

In order to estimate the sum in the r.h.s. of (33), we first note that coherent contributions come 
only from atoms whose z-position satisfies the condition kozj/i2r) < 1. The coherent terms 
appear at the Bragg angles 



sm\6g/2)^^, q^\,2,.:qmax (34) 
koa 



Expanding the terms in the phase of Se/f depending on 6 around one Bragg angle we find 

Seffiko,r,t)~iD/a) ^ gao[sine,(e-e,)«.«+4^] (35) 

where n^^ ~ [(l/a)(^)'^^]. Exploiting the foil's symmetry in the z direction we approximate the 
sum in (35) as 

due'-^--\ due'^kldn^0ci{e-0cifu^] (36) 
2 Jo 

where Umax ~ (r/ko)^^^. An explicit formula for the above integral can be given in terms of 
error functions. However, a qualitative imderstanding of its behavior is ofl'ered by the approx- 
imation ^ 1 + whereby it follows that due'^ ^ e'^^^^{Umax - I'-yuLix) 
jj due^!~u^^ = ^u^^e^^ . Substituting the above expressions in (36) we find 



Seff{ko,r,t) ~ 2(£)/a)(l/a)e^[M„„,- i^M^^- i^2^in2 0,(0-0,)2„3^J (37) 



Taking into account also the difl'use term, the final form of the outgoing wavefunction is 

^cu,,oin,- 2?^^^^iD/a)e-'^fir,0)e''»^ 
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Figure 6: (a) Local deformation of tfie separator at tlie time t = Iq/vq after tlie inclusion of the Bragg angles 6^ (Eq.34) in 
the effective Fraunhofer function (we mark the first three angles as A, B, C). (b) The deflection of the Bohmian trajectories 
at the channels of radial flow formed around the Bragg angles A,B corresponding to di = 0.23... and O2 = 0.33... 
respectively), (c) The concentration of the scattered trajectories close to Bragg angles in a larger scale, (d) Angular 
distribution corresponding to the numerical trajectories of (c). The dashed lines denote the exact positions of the Bragg 
angles. 
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where the sum is considered with respect to all Bragg angles, while the following estimates hold 
for the functions Uq and O^: 



(39) 



and 



<bq ~ tan 



-1 



^-3 + irfeosin^ 61^(0 -0^)2 



) 



(40) 



sufficiently far from the target. The last equation implies that at angular distances \6 - 6q\ ~ 
7T/(rk()y^^, the particles' Bohmian trajectories acquire a transverse velocity v, - (1 /r)d^q/d6 
pointing towards the direction of the straight line with inclination equal to tan 8q, while v, = 
exactly at d = 6q. Furthermore, the presence of the coherent terms in Seff causes a local 
deformation of the separator around the Bragg angles, as shown in Fig. 6a. We note that the 
separator comes locally closer to the center at the directions corresponding to the Bragg angles, 
since the magnitude of iff outgoing is locally enhanced due to the local peaks of the functions Uq. 

The effect of this deformation on the Bohmian trajectories is analogous to the one described in 
[51]. Namely, this deformation results in the formation of local channels of radial flow, whereby 
the Bohmian trajectories are preferentially scattered around the Bragg angles. An example of 
this concentration is shown in Fig. 6b. Clearly, the inclusion of the coherent terms causes a 
variation of the angular distribution of the Bohmian trajectories, by creating local maxima of 
the density around the Bragg angles (Oi = 0.23..., 62 = 0.33... in Fig.6b). Figure 6c shows 
this concentration in a larger scale, while Fig.6d shows the angular distribution corresponding 
to the trajectories of Fig.6c. This distribution exhibit clear peaks at all the angles - 0q (the 
first local maximum around = 0.1 is not due to a concentration at a Bragg angle, but it is 
only caused by the trajectories moving nearly horizontally, i.e. within the support of the ingoing 
wavepacket). We note that plots of the quantum trajectories in a different scattering problem 
(atom surface scattering), appearing in [21][22], show a similar qualitative picture as in Fig.6b, a 
fact which was identified in that case too as a dynamical effect of the quantum vortices. We may 
thus conjecture that the quantum vortices play an important role in a wide context of different 
quantum-mechanical diffraction problems. 

Finally, it should be stressed that the modification of the outgoing wavefunction according to 
Eq.(38) only influences the Bohmian velocity field in the transverse direction, while the radial 
flow of all Bohmian trajectories (as in Fig.6) takes place at a constant speed hkajm. Thus, the 
emergence of a diffraction pattern does not influence estimates on times of arrival or the times of 
flight of the particles to detectors placed at the same distance from the center, independently of 
the angle 6. This subject is now discussed in section 4. 

4. Arrival times and times of iliglit 

An important practical utility of the quantum trajectory approach regards the possibility to 
unambiguously determine the probabihty distributions of the so-called arrival times, or of the 
times of flight of the scattered particles. This question is of particular interest, because it is related 
to the well known 'problem of time' in quantum theory (see [55, 56] for reviews). This problem 
stems from a theorem of Pauh [54], according to which it is not possible to properly define a 
self-adjoint time operator consistent with all axioms of quantum mechanics. This implies that 
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the usual (Copenhagen) formahsm based on state vectors or density matrices is not appHcable 
to a quantum-theoretical calculation of probabilities related to time observables. In fact, in both 
Schrodinger's and Heisenberg's pictures, time is considered only as a parameter of the quantum 
equations of motion. 

Among various proposals in the literature aiming to remedy this gap of standard quantum 
theory ([55]), the Bohmian formahsm offers a straightforward solution. This is in principle sub- 
ject to experimental testing, as will be proposed below. Furthermore, as was mentioned in the 
introduction, the use of a wavepacket approach allows for a comparison of the Bohmian approach 
with two other main approaches to the same subject, namely the 'history approach' (based on 
Feynman paths) [57] [58], and the Kijowski approach, based on so-called 'Bohm-Aharonov oper- 
ators [59]. We should emphasize, however, that so far in the literature the latter approaches were 
given a consistent formulation only in the case of asymptotically free wavepacket motion [55], 
while their implementation in the case of scattered wavepackets is an open issue. This will be 
discussed in a future work. In the present paper, on the other hand, we focus on results regarding 
the time observables as defined in the Bohmian approach, and only provide a rough comparison 
of what should be expected by other theories of time observables. 

Any definition of a time observable in quantum mechanics (and also in classical physics) 
requires the occurrence of two events serving as the 'start' and the 'stop' event in the process of 
timing. 

In the definition of the arrival times of particles to detectors, we take as start event the prepa- 
ration of the whole initial wavefunction at a certain moment f = f i , which can be conveniently 
set equal to t\ - 0. The stop event is the detection of the particle at a time t = t2 > h. The arrival 
time of the particle to the detector is tamvai - h- t\. 

The above definition of arrival times is independent of the adopted picture of quantum me- 
chanics, since its only requirement is to assume that the preparation of an initial quantum state 
can be controllable in time (see [71]), i.e. that replicas of the same state can be prepared at 
any given time t\ (a technique realizing this experimentally will be proposed below). However, 
a consistent calculation of the arrival time probabilities by the various pictures is still an open 
theoretical issue (see [55]). 

On the other hand, the definition of a time of flight for a particle depends on the adopted 
picture of quantum mechanics, since it requires the use of some notion of spacetime paths that 
the particles presumably follow within the picture's framework. In fact, the time of flight is 
defined as the time elapsing between the crossing by the particle of two surfaces S i and 5 2 in the 
configuration space. Thus, this time is different from the arrival time, and the difference depends 
on where exactly the particle lies within the support of the initial wavefunction. 

In the case of particle diffraction, it is convenient to choose the two surfaces as shown in 
Fig.l: S\ is taken normal to the z-axis at a point z = -lo, while 5 2 is a spherical surface sur- 
rounding the target at the radial distance r2-h. The time of flight T{6) depends on the scattering 
angle 6, and the quantity of interest is the difference T{02) - T{0\) for two angles 9\, 62- This 
difference is independent of /q, provided that /q is such that both 5 1 and 5 2 are sufficiently far 
from the target. 

We now discuss separately the arrival time and the time-of-flight probabihties in the setup of 
Fig.l. 

4.1. Arrival times 

The motion of all scattered particles (like in Fig.4) in the domain beyond a sphere of radius 
2/0 around the center can be considered as a radiaUy outward motion with constant speed, since, 
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by taking ^ 4'oiiigoing (where 4'outgoing is given by Eq.(16)), the Bohmian equations of motion 
in spherical coordinates read dridt = vq = hko/m, dOjdt - 0. The last equation is modified 
when the diffraction terms in Seff are taken into account. This modification, however, does not 
influence the motions in the radial direction, which, as shown now, are the only ones afifecting 
the distribution of arrival times to a detector placed at a distance Id and at any fixed angle 6. By 
definition, the latter distribution is given by 

P.rr.am = (41) 

where ANgji^(t) denotes the number of particles within the (assumed fixed) detector conic aper- 
ture dD.D = sin 6A9r)A(/)r> around the angle arriving to the detector between the times t and 
t + At. Since in the Bohmian approach all the particles move with constant speed, we have 

^Ngj^jt) ANej^(t) Ar . . .UKq 

At Ar At m 

In the case l» D » a, substituting Eq.(22) and making the usual approximations » ht/m, 
» ht/m results in 

-(Id*Io-W 

Parrivalit) ^ PqC ^ (42) 

where Pq is a normalization constant. On the other hand, in the case D » I » a (via Eq.(25) 
we find 

ParrivaM = P'qC »in^»^ (43) 

The main result can be summarized as follows: in either case l» D or D » I, the arrival time 
distribution is a localized distribution (Gaussian, around the mean time (£d + ^o) / vq), whose dis- 
persion is always of the order of Vq ^ x the maximum of the transverse and longitudinal coherence 
lengths. 

The latter property implies that the trajectory approach makes predictions regarding the ar- 
rival time distribution which depend on two main beam parameters, and are thus testable in 
principle by concrete experimental setups. One possible proposal in this direction is the use of 
the so-called laser-induced cold field emission technique (see [72]). In this technique, a cold- 
field electron source (nanotip) is exposed to well separated in time focused weak laser pulses of 
time width ~ 10-lOOfs. The photo-emitted electrons are accelerated towards the anode, whereby 
their initial state can be effectively described by an ingoing wavefunction of the form (15). The 
scattered particles pass through detector placed at fixed angles 6 as indicated in Fig.l. The key 
point to notice is that time measurements in such a setup conform with the definition of the ar- 
rival times, since a detection of the triggering laser pulse can serve as a start event marking the 
initial time when the whole state was prepared, while a later detection of a scattered particle 
serves as the stop time t2. We propose that by monitoring the electron beam one can achieve 
different values of / and D, thus probing quantitatively the predictions for Pamvaiit, 6) as given 
by the quantum trajectory approach. 



4.2. Tlmes-of-flight 

The total time of flight from a point on 5 1 to a point on 5 2 is a function of the initial conditions 
iZQyRo)- Using the information from the swarm of numerical Bohmian trajectories of Fig.4, this 
function can also be quadratically interpolated by the numerical data on grid points. The mean 
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time T{6) for all initial conditions leading to the same 6 is then found numerically. Since the 
choice of surfaces S \ and 5 2 in Fig. 1 is arbitrary, the invariant quantity of interest is the difference 
T{6) - T(do), where 60 is a fixed reference angle. Figure 7 shows this difference for the numerical 
trajectories of Fig. 4a. 

To estimate T(0) theoretically, we make use of Eq.(32), yielding the locus £.{6) of all initial 
conditions leading to a scattering close to the angle 6. Using also the separator equation (Eq.27), 
we also find the point (zs,Rs) where the moving separator encounters an orbit moving horizon- 
tally from (z£(g){R),R) at f = 0, with speed vq. Thus, we set R = R, and Zs = Rs/ tan(0). The 
time of flight of this trajectory from 5 1 to52 is then t(zj:(e)(R),R) = (Zs + 2/o - h - R/ smd)/vo. 
The mean time of flight T(d) can then be approximated by T{6) » JT,^^^ 2nR\tffi„{z£(e)iR),R, t = 

0)\\zj:(,g)iR),R)dR. We flius find: 

AT = 7(61) -7(62) » ^[taii(02/2))-taii(0i/2)] (44) 

Vo 

whereto = [ V21n(Co) + 1/(1 + V21n(Co))], wifli Co = 8neoklh^/i\ZiZ\q^^mp^/^d^'^). 




Figure 7: The time difference T{9) - r(150°) (see text) for the Bohmian trajectories of Fig.4. The smooth soUd curve is 
the theoretical prediction of Eq.(44) while the dots represent numerical results. 

Due to Eq.(44), the time difference T(9i) - T{02) has an 0(D tan(0/ 2) /vq) dependence on 
0. In fact, it is noticeable that Eq.(44), which gives the difference of the mean times of flight in 
the wavepacket approach, turns to be identical to the estimate of [51] (their Eq.(26)), referring 
to the plane wave approximation. This is expected, since the plane wave approximation can be 
considered as a Umiting case of the wavepacket approach with / » D, corresponding to the limit 
/ 00. 

One more interesting remark concerning the times of flight found by the Bohmian approach 
is that the difference T(0i) - T(02) predicted in Eq.(44) has a completely different behavior from 
analogous quantities calculated in the framework of other theories of quantum time observables. 
While we defer a detailed reference to this problem to a future work, here we give some rough 
estimates concerning the sum-over-histories and the Kijowski approaches referred to in the in- 
troduction. We find 



AT = TiOi) - T{e2) « ^ In 

(in sum-over-histories formalization, semiclassical approximation) 
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l l +cot^(6li/2) 



(45) 



and 

AT = T(6i) - r(6»2) = (in Kijowski formalization) (46) 

An outline of the derivation of these formulae is given in Appendix I. 

A comparison of all three approaches yields that (i) the sum-over-histories approach (which, 
using a semiclassical approximation yields essentially the same result as in classical scattering 
theory), predicts a mean time difference depending on the particle velocity vq, by the scaling law 
AT ~ [me/m][10'^m sec^Vvo]^ ■ lO^^^^ec. (ii) The Kijowski formalism predicts no time differ- 
ence, (iii) The Bohmian formalism predicts a time difference depending on vo as well as on the 
transverse quantum coherence length Z). The scaling AT" ~ [Z)/lyumJ[10**m sec^'/i'ollO^'^sec 
holds. 

As a final conclusion, we propose that experiments aiming to measure time observables in 
setups of particle diffraction may provide new insight into fundamental problems such as the role 
of time in quantum mechanics. In particular, the predictions of the de Broglie Bohm theory are 
within the possibilities of present day experimental techniques. 



5. Semiclassical limit (Rutherford scattering) 

So far we have considered charged particles with quantum coherence lengths much larger 
than the distance between nearest neighbors in the target. However, this study does not cover 
the so-called short wavelength limit, as e.g. in the case of a-particle or ion scattering. In this 
case, the quantum coherence length becomes comparable to or smaller than the distance between 
nearest neighbors in the target. As a result, such particles 'see' each of the atoms in the target 
as an individual scattering center and they do not interact with the target lattice as a whole. 
Furthermore, both D and I become comparable to a classical 'impact parameter' b (of the order 
of a few fermi) which is relevant to the classical description of Rutherford scattering. 

The incorporation of h in the wavefunction model can be done essentially as described in 
[69], assuming a Gaussian form of his wavepacket, referred to by the notation and aligning his 
vector denoted by b along the x-axis of our coordinate system we have: 

^) = 7^ r + A(0))— (47) 
KLnyi^ J r 

with 



c(k) = -^.<"-i'— 1— . e-'"'" (48) 
and 

Z\Zqi m 1 

fm = - -2 ^ . . . „ (49) 

After the calculation of the above Gaussian integrals we are led to the following wavefunction 
model: 

if/ir, t) = (//ingoingir, t) + ll/outgoingir, t) (50) 

where 



l^mgoingiVJ) = Aexp - 
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4', 



outgoing 



(r,0 



I 47reo 



exp 



2ffkl sm\0/2)r^ 



(52) 



2(D2 + ihtim) 



+ Hkor - hk^t/lm) 



where 



A = 



D 



1 



fl/2 
)tT77(- 



1/2 



(53) 



is a nearly constant quantity, not affecting the Bohmian trajectories, and vq = hko/m. The time 
f = in the above formulae is taken so that, in the absence of scattering, the center of the ingoing 
wavepacket crosses the plane z = at the moment f = 0. Furthermore, we consider the Bohmian 
trajectories at positive or negative times satisfying \t\ < mD^/h, i.e. smaller than the decoherence 
time of the packet. 

The outgoing term is modulated by the Gaussian factor 

(r - \'o/)- + h- 



^^^\ 2(D^ + iht/m) 

This factor implies that a replica of the ingoing wavepacket propagates from the center outwards 
as a spherical wavefront of the outgoing wave, albeit by a phase difference ivot with respect 
to the ingoing wavepacket. This new factor is the most important for the analysis of Bohmian 
trajectories, because it implies that the form of the latter depends crucially on the choice of the 
value of the parameter b, which actually changes the form of the wavefunction. 




Figure 8: Three bohmian trajectories guided by the wavefunction defined in Eqs.(51) and (52), with D = 1 = lOfermi, 
Zi = 2, m = 7.1 X lO^mf (alpha particle), Z = 79 (gold), ko = 10^'^m"', and three difterent impact parameters 
bi = lOfm, &2 = 12fm, = 15fm. All three trajectories are initially posed at the centers of their corresponding guiding 
wavepackets. 



A careful inspection of Eqs.(5 1) and (52) shows that the spherical wavefront emanating from 
r = encounters the ingoing wavepacket at a time tc which decreases as b decreases. As a 
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result, the Bohmian trajectories, which are forced to follow the motion of the radial wavefront 
after the collision, are scattered to angles which are larger on the average for smaller b. Thus, the 
Bohmian trajectories recover on the average the behavior of the classical Rutherford trajectories. 
An example of this behavior is given in Fig. 8, showing three Bohmian trajectories corresponding 
to an initial condition taken at the center of the wavepacket defined by Eq.(51) at the time / = 
-D^m/H, ensuring that the packet's spreading does not change appreciably in time up to the 
moment of collision with the outgoing wavefront. The three Bohmian trajectories cross each 
other, yielding a larger scattering angle for a smaller initial distance from the z-axis, i.e. they are 
close to the familiar classical picture. It should be noted, however, that this closeness is only in 
an average sense, since the exact form of a Bohmian trajectory guided by a wavepacket depends 
on where exactly the initial condition of the trajectory lies with respect to the center of mass of 
the initial packet. In fact, for one fixed value of b one obtains a swarm of de Broglie - Bohm 
trajectories (with initial conditions around this value of b). These trajectories are scattered in 
various directions and they do not cross each other, while they can define (in a statistical sense) 
a most probable scattering angle. However, this angle increases as b decreases. Hence, we 
conclude that when we consider the 'semiclassical limit' of small wavelengths as well as small 
quantum coherence lengths, the Bohmian trajectories yield results which agree on the average 
with the classical theory of Rutherford scattering. 

6. Conclusions 

We developed a wavefuntion model providing a wavepacket approach to the phenomenon of 
charged particle diffraction from thin material targets, and we employed the method of the de 
Broglie - Bohm quantum trajectories in order to interpret the emergence of diffraction patterns as 
well as to calculate arrival time probabilities for scattered particles detected at various scattering 
angles 6. Our main conclusions are the following: 

1) In both cases when the longitudinal wavepacket coherence length / is larger than the trans- 
verse wavepacket coherence length D, or vice versa, the outgoing wavepacket has the form of a 
pulse propagating outwards in all possible radial directions, with a dispersion cTr which is of the 
order of the maximum of / and D. Furthermore, in the case D » I (applying e.g. to cold-field 
emitted electrons), ar depends on as o"r ~ sin 6D. 

2) We study the structure of the quantum currents in the above model. We provide theoretical 
estimates regarding the form and time evolution of a locus called separator, i.e. the border 
between the domains of prevalence of the ingoing and ourgoing quantum flow. We show how the 
separator forms channels of radial flow close to every Bragg angle, and leading to a concentration 
of the quantum trajectories to particular directions giving rise to a diffraction pattern. 

3) The deflection of quantum trajectories is due to their interaction with an array of quantum 
vortices formed around a large number of nodal points located on the separator. We show exam- 
ples of the quantum flow structure forming a 'nodal point - X-point complex' around any nodal 
point, and we calculate the form of the stable and unstable manifolds yielding the local direc- 
tions of approach to or recession from an X-point. In view of the similar role played by quantum 
vortices in different examples of difiraction problems [21] [22], it can be anticipated that the 
mechanism of emergence of the diffraction pattern described in section 3 is quite general. 

4) We compute arrival time probability distributions for both cases l» D and D » I using 
the de Broglie - Bohm trajectories of particles detected at a fixed distance and various scattering 
angles with respect to the target. In all cases, the dispersion of the arrival time distribution turns 
to be cr, ~ Vq ' X max(/, D), where vq is the mean particles' velocity. We propose a realistic 
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experimental setup aiming to test this prediction for electrons. We also calculate time-of-flight 
differences, where the time of flight is defined as the time interval separating the crossing by a 
de Broglie - Bohm particle of two fixed surfaces located in the directions of asymptotically free 
motion before and after the target. We discuss the ambiguity of the definition of the times of flight 
when using different approaches besides de Broglie - Bohm, and provide a rough calculation in 
the framework of the sum-over-histories approach and the Kijowski approach. 

5) We finally examine how the de Broglie - Bohm trajectories recover (in a statistical sense) 
the semiclassical limit of Rutherford scattering, by examining the form of the quantum trajecto- 
ries when the packet mean wavelength, as well as / and D become smaller than the inter-atomic 
distance in the target. In particular, we incorporate an impact parameter b in the wavefunction 
model and demonstrate that the de Broglie Bohm trajectories are scattered on the average at 
larger angles 6&sb decreases. 
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Appendix A. Time-of-Flight differences outside the Bohmian formalism 

We hereby outline the derivation of Eqs. (45) and (46) estimating the time-of-flight differ- 
ences T(92) - T{9i) by the sum-over-histories and the Kijowski approach respectively. 

i) Sum-over-histories approach. Neglecting the question of consistency (see [55]), a rough 
calculation of the difference T(02) - T{9\ )) in the sum-over-histories approach can be made in the 
framework of the method developed in [58J. Denoting by Q a fixed space-time volume with time 
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support within the interval from two fixed times ?a to ts, the probabiUty of a particle crossing Q 
is 



dhA<^(B;n;A)ilf(rA,tA) 



(A.1) 



where A = (rA,tA),B = {ybJb), while 0(B;Q;A) = SysfB^n-A) e''^'^''^*^^ is the sum over all 
Feynman paths y with fixed ends A and B passing through Q. The quantity P{tA, tg, r^; f2) - 
I ^ dvA^iB; Q; A)ij/{rA, Ia) P ^^^r^ is identified as the probability that an electron being anywhere 
in space att = Ia reaches a volume + (Ptb att = ts by passing first through Q. Let Si be 
a fixed surface normal to the beam's central axis (Fig.l) at -Zq < Zs, < 0. We choose Q by the 
conditions that belongs to an area element A5 1 on 5 1 , around the point R = 0,zq = Zsi, while 
to < tn < to + At, where to > tA- Let now 5 2 be a spherical surface or radius /q around O, and 
a point on 52 in the plane of Fig.l. Setting d^rg = A^avoAf, where AS 2 is an area element 
on 5 2 around r^, the mean time of flight from AS 1 to A52 becomes a function of 6 only, given 

by T(0) = / (tB - to)P{tA = 0, tB, rsiO), n)dto where Po = [J PHa = 0, tB, r^Ce), Q)Jfo)"\ To 
estimate P(tA, tg, r^; Q.) we extend results of Hartle [57] and Yamada and Takagi [58] in our case. 
For (^(B; £1; A) we adopt a 3D extension of a formula proposed in [58]: 



<D(B;Q;A) 



-HI 



dh^{rB,tB\r',to + At) 



0(r', to + At\ Q; r, to)^{Y, to; ta, tA) 



(A.2) 



where <I)(r, fo; r^, f^) is approximated by a free Feynman propagator, 0(rB, f^; r', fo + At) = 

J d^ke 2^ (pk(x')0k(XB) (with 0k as in Eq.(ll)), while (^(r',to + At;Q.;r,to), is approxi- 
mated by a 3D analog of Hartle's approach [57] 



<l)(r',fo + Af;n;r, fo) 



Wo+Af r ^ / ^ \3/2 

i '^\L.'^'^''[2nihito^At-t)l 



•Si I \ R|^2(fo+Al-0^ 2(Mo) 



■to 



(A3) 



where r^j are points on A5i and Rsi = - zsi^z- An exact calculation of all integrals is 
untractable. Through a stationary phase approximation, however, one obtains for fixed < 
and if/ir, tA) - ^'m, that P(tA, ts, r^; Q.{to)) is peaked essentially around a mean 'classical' time of 
flight for a trajectory starting from the center of the wavepacket |^,„|, scattered by an atom at O, 
and arriving to a point on AS 2. For two angles 6\, 62, the mean time of flight difference can be 
estimated as 



T(0i) - T(62) 



In 



jl +COt2(0i/2)' 
V 1 + COt2(6'2/2) 



(A.4) 



i.e. we find Eq.(45). 
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ii) Kijowski approach. The Kijowski approach [59], assuming one-dimensional wave-packet 
propagation along some direction z, we set 




1/2 



2 



c{k)e' 



'■T/2m+ikz 



(A.5) 



to be the probability that the particle arrives on a normal surface at a point z between times T 
and T +dT (c(k) is the wavefunction in momentum space). If c(fc) 

^^-(k-kof/2^l-ikzo is a narrow 
packet (^0 » cy); neglecting exponentially small negative component terms of cik), we find 



i.e. 11(1, z) practically coincides with the flux function J{z,T) = (hko/m)\tJ/ix,T)\^. We apply 
Kijowski's formaUsm in the setup of Fig.l for the asymptotically free wavepacket motions at 
times long before or after t = Iq/vq. The mean time-of-flight from 5 1 to a point of fixed on 5 2 
can be written as T(0) -<t2> - <ti >, where t\,t2 are the arrival times to S 1 and 5 2. Applying 
(A.6) we find < t\ >= Ii/vq, < ti >= 2/o/vo. Thus, T(9) = (2/o - /i)/vo independently of the 
scattering direction, i.e. 



UiT,z) = ^ l+0((r||Ao) \<ffiz,T)\' 



(A.6) 



noi) - 7(02) = 



(A.7) 



i.e. we find Eq.(46). 



27 



